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ABSTRACT 

We present recent improvements in polarization analysis with the INTEGRAL SPI data. The SPI 
detector plane consists of 19 independent Ge crystals and can operate as a polarimeter. The anisotropy 
characteristics of Compton diffusions can provide information on the polarization parameters of the 
incident flux. By including the physics of the polarized Compton process in the instrument simulation, 
we are able to determine the instrument response for a linearly polarized emission at any position angle. 
We compare the observed data with the simulation sets by a minimum \ 2 technique to determine the 
polarization parameters of the source (angle and fraction). We have tested our analysis procedure 
with Crab nebula observations and find a position angle similar to those previously reported in the 
literature, with a comfortable significance. Since the instrument response depends on the incident 
angle, each exposure in the SPI data requires its own set of simulations, calculated for 18 polarization 
angles (from 0° to 170° in steps of 10°) and unpolarized emission. The analysis of a large amount of 
observations for a given source, required to obtain statistically significant results, represents a large 
amount of computing time, but it is the only way to access this complementary information in the hard 
X-ray regime. Indeed, major scientific advances are expected from such studies since the observational 
results will help to discriminate between the different models proposed for the high energy emission 
of compact objects like X-ray binaries and active galactic nuclei or gamma-ray bursts. 
Subject headings: gamma rays: general - instrumentation: polarimctcrs - methods: data analysis - 
polarization 



1. INTRODUCTION 

Throughout the universe we observe powerful engines 
which accelerate particles to immense energies. The pre- 
cise details of how these engines function are still poorly 
understood, but polarization measurements of the high 
energy radiation are pivotal to gain a deeper insight into 
the physical environment in these systems. 
Highly ordered geometries are generally required to pro- 
duce a net polarization signal in the emission from a 
source, and strong magnetic fields are usually thought 
to play a leading role in the emission processes. The 
strength and direction of a polarization signal will lead 
back to an understanding of the physical mechanisms 
and their geometry at the emission site. There have 
been few attempts at measuring hard X-ray polarization 
since measurements are hampered by large backgrounds 
and systematic effects within the detectors. However, 
the advent of modern fast computing clusters has made 
large scale simulation of an instrument's response now 
possible. By combining instrument data and results 
from detailed Monte-Carlo Mass-Model simulations us- 
ing GEANT4, it is possible to put constraints on the 
polarization characteristics of the hard X-ray flux emit- 
ted by a source. Using the Compton scattered events 

* BASED ON OBSERVATIONS WITH INTEGRAL, AN ESA 
PROJECT WITH INSTRUMENTS AND SCIENCE DATA CEN- 
TER FUNDED BY ESA MEMBER STATES (ESPECIALLY 
THE PI COUNTRIES: DENMARK, FRANCE, GERMANY, 
ITALY, SPAIN, AND SWITZERLAND), CZECH REPUBLIC, 
AND POLAND WITH PARTICIPATION OF RUSSIA AND USA. 

1 Universite de Toulouse, UPS-OMP, IRAP, Toulouse, France; 

2 CNRS, IRAP, 9 Av. Colonel Roche, BP 44346, F-31028 
Toulouse Cedex 4, France 

3 Createc Ltd, Unit 8, Derwent Mill Commercial Park, Cocker- 
mouth, Cumbria CA13 0HT, UK 



in the SPI instrument on INTEGRAL (130-8000 keV), 
constraints have been already put on the percentage po- 
larization in t he prompt gamma-ray flux of gamma-ray 
bur sts ( GRBs: iMcGlvnn et al.ll20"07j ) and the Crab pul- 
sar (|Dean et al.ll2008l ). The latter showing a remarkable 
alignment to the rotational axis of the spinning Neutron 
star. 

However, improvements on the methods and techniques 
employed in the search for polarization could be achieved. 
This paper discusses new developments that have been 
made in the mass-model and in the analysis procedure 
for the SPI instrument. 

2. DETECTION OF POLARIZED EMISSION 

The signature of polarization in the source signal can 
be found through the interaction process of photons 
within detectors. 

In the hard X-ray regime, the dominant process is Comp- 
ton scattering. In the case of a linearly polarized flux, the 
azimuthal angle (j> of the scattered photon is no longer 
isotropic. The probability for the photon to be scattered 
at a polar angle 9 and an azimuthal angle <f> is given by 
the Klein-Nishina differential cross-section 
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where r$ is the classical electron radius, e is the ratio 
between the scattered and incident photon energies, and 
4> is the azimuthal scatter angle defined as the angle to 
the polarization unit vector. An azimuthal anisotropy is 
observed in the scattering direction, producing a greater 
number of scatters in a direction perpendicular to the 
emission polarization angle (P. A.; see Figure [T]). This 
formula shows that as the energy increases the angu- 
lar distribution related to the polarization becomes more 
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Figure 1. Klcin-Nishina differential cross-section for a photon 
polarized parallel to the x-axis, seen from the direction of incidence, 
showing the anisotropy in the azimuthal scatter direction. As the 
energy increases the photon is more likely to be forward s cattered 
and the distribution becomes more isotropic (Clark 2007). 

isotropic. 

Using a suitable arrangement of detector elements, this 
anisotropy can be used to determine the direction and 
degree of polarization of any observed emission. 

3. SPI AS A POLARIMETER 

SPI is a coded-aperture telescope on board the IN- 
TEGRAL observatory operating between 20 keV and 8 
MeV. The description of t he instrument and its per- 
formances ca n be fo und in iVedrenne et al.l (|2003[ ) and 
iRoaues et all (|2003[ ). SPI consists of 19 hexagonal Ger- 
manium detectors measuring 5.6 cm flat-to-flat and 7 
cm deep, tessellated into a hexagonal shape (see Figure 
[2|). The center-to-center distance between detectors is 6 
cm. Originally designed for fine spectroscopy, SPI can 
be used as a Compton polarimeter. The measurement of 
polarization using Compton events relies on the positions 
of the initial and final photon interaction in the detec- 
tors. In SPI, there is no positional information available 
inside a detector so only the multiple events where the 
photon deposits its energy in more than one detector can 
be used. This kind of event is identified by energy de- 
posits occurring within a coincidence window of 350 ns 
and not detected by the anti-coincidence system. The 
effective ar e a of SPI for multiple events is reported in 
lAttie et all (|2003l ) and range from 10 cm 2 to 50 cm 2 be- 
tween 150 keV and 1 MeV. Since INTEGRALS launch, 
four detectors have failed: detector 2 (2003), detector 17 
(2004), detector 5 (2009), and detector 1 (2010). These 
failures have resulted in a decrease of the effective area 
of the instrument down to ~ 52% of the original area for 
multiple events. 

Besides its effective area, the quality of a polarimeter de- 
pends on the sampling of the Compton scattered events 
distribution (Figured]) by the detectors. An usual figure 
of merit for polarimeters is the Q-factor, defined as the 
ratio of the parallel and perpendicular components of the 




Figure 2. The 19 germanium detectors before installation into 
the spacecraft (Credit: ESA). 

scattered distribution. The Q-factor, which is energy de- 
pendent, has been calculated fo r SPI to be ~ 0.24 over 
the energy range of 0.1-1 MeV (jKalemci et al.l [200 71) . 
When searching for polarization in a source emission, any 
scattered distribution on the detector plane will be mod- 
ulated by the coded mask shadow and the dead detectors. 
These effects can be mistaken for a polarization signal if 
they are not properly accounted for. Instead of looking 
for an anisotropy in the scattered distribution, a much 
simpler method is to simulate the detector response for 
a polarized flux and compare it to the measured data. 
In this way, polarization analysis with SPI is perform 
in the same way than imaging but using the instrument 
response with one additional dimension: the P.A. Thus, 
for each SPI observation we need simulations for different 
P.A.s in the source emission. 

4. SIMULATING THE SPI RESPONSE FOR A POLARIZED 
EMISSION 

4.1. The Geant4 Model 

The presented simulations are based on a model re- 
coded f rom the GEANT3 TI MM (for The Integral Mass 
Model; iFerguson et~aT1 12003D into the GEANT4 fr ame- 
work. GEANT4 (|Geant4 Collaboration et"al|[200l uses 
object orientated C++ language and includes extra 
physics such as polarization and material activation. The 
Geant4 software package has been tested a gainst calibra- 
tion d ata from a prototype detector array (jMizuno et al.l 
2005) and reliably reproducing the instrument response 
to a polarized beam. In the new model of INTEGRAL, 
the SPI and JEM-X detectors have been both fully mod- 
eled. However, at the current time, the IBIS detector is 
only described with an approximate geometry (Figure 
This is not a problem for our polarimetry studies, since 
the work revolves around the SPI instrument. Masses 
situated away from the SPI instrument will have little or 
no impact on the results. Conversely, the SPI instrument 
has been modeled with a high degree of accuracy. It com- 
prises the active veto, the Ge detector elements inside Al 
capsules, surrounded by the Be housing, and the mask 
assembly. The GEANT4 incarnation of the INTEGRAL 
model was specifically written to assess the background 
in the SPI instrument. This means that several compo- 
nents of the geometry and particle generation processes 
have not been considered but could become important 
when looking at the emission from an astronomical ob- 
ject. For example, the exact geometry of the coded mask 
and the ability to generate photons with a given spectral 



Hard X-ray polarimetry with INTEGRAL SPI 



3 



shape had to be implemented. 

A great deal of effort has been spent on updating the sim- 
ulation to produce a much more accurate model of the 
SPI instrument, including several missing components in 
the geometry. In the current status, our simulations take 
into account the effects of non-functioning (dead) detec- 
tors, off-axis illumination and mask shadowing . Com- 
pared to the model version used by iDean et al.l (|2008j ). 
we have included the degraded transparency of the cen- 
tral mask pixel which is 66% at 130 keV and 80% at 
600 keV. We have also modified the anti-coincidence sys- 
tem. In fact, the most important part surrounding the 
detectors was not active in the code and we introduced 
the anti-coincidence thresholds (^100 keV). These mod- 
ifications have an important impact on the number of 
events in the outer detectors and drastically change the 
instrument response. The individual germanium crys- 
tal dimensions have been slightly reconsidered (result- 
ing in a few % smaller volumes) modifying the instru- 
ment sensitivity. In particular, the misalignment of the 
SPI telescope axis relatively to the INTEGRAL satellite 
pointing axis has been corrected, which corresponds to 
a systematic of ~ 0°.15 in the pointing accuracy. This 
correction is essential for the simulation of the instru- 
ment response because it modifies the mask shadowing 
and thus, the event distribution over the detectors (~ 6% 
of the events). These important changes in the instru- 
ment simulation lead to a good agreement with the SPI 
standard response (see Section |4~4|) . 
The polarimetry capability of the SPI instrument is 
based on the Compton efficiency of the detector plane, 
i.e., the non-diagonal terms of the response matrix (~ 
20% of the total efficiency) . This Geant4 model allow us 
to carefully evaluate them, through their dependence on 
the P.A. of the incident radiation. It provides a robust 
technique, using the whole instrument response and thus 
giving the best sensitivity to the polarization character- 
istics. 

4.2. Simulated Data Production 

The simulated source is characterized by its direction 
in equatorial coordinates (R.A., deck), an energy spec- 
trum, and a P.A. This angle corresponds to the position 
angle of the electric vector measured from north to east 
on the sky. For each pointing, simulations are carried 
out at 18 PA.s (0°-170° in 10° steps) and 1 unpolarized. 
Only the P.A.s between 0° and 170° need be modeled due 
to the 180° symmetry of the Compton scattering process. 
To properly simulate the photons distribution, the data 
is pre-analyzed and fitted with a simple spectral shape. 
This is performed with the c lassical SPI sp ectral analysis 
using the XSPEC software (|Arnaudl ll996fi. The spectral 
models form the basis of a probability distribution for 
the fired photons so that a single photon may be sim- 
ulated at any energy, but after millions of photons, the 
spectrum should resemble that of the observed data. In 
theory, it should be possible to simulate a flat input spec- 
trum with a fixed number of photons in each energy bin 
and to normalize a posteriori the number of photons in 
each energy bin. Simulations would need to be carried 
out only once and spectral changes could be handled in 
the analysis. However, this will needlessly simulate many 
photons at high energies and would drastically lengthen 
the running time of the simulations. 




Figu re 3. The GEANT4 INTEGRAL model geometry IClark I 

Besides, a small difference in the input spectrum will not 
lead to a big effect in the simulation. We have compared 
two simulations, with a power law slope of -3.0 and -1.5, 
respectively. This results in differences less than < 10% 
in the count rate distribution recorded in the detectors. 
For each simulation, millions of photons are fired from 
the direction of the source. Their initial positions are 
randomly distributed on a 140 x 140 cm surface, 500 cm 
away from the detector plane, so as to illuminate all the 
instrument. This area has been chosen large enough to 
include the direct lighting from the source and the scat- 
tered events from the structure in the detected events. 
For a typical power law spectrum with a slope of -2.2, 
50 million of photons are processed for each simulation 
resulting in ~ 500,000 single events and ~ 50,000 mul- 
tiple events in the detector. For a source flux of 1 Crab 
observed during one pointing (~2 ks), this provides an 
oversampling of approximately 40 times. 
For one photon fired (event), all the energy deposits in- 
side the detectors are stored if no deposit occurred in the 
anti-coincidence system. 

The number of simulations needed makes this part of the 
analysis very time consuming. One simulation takes ~ 6 
hr to run on a single processor. To produce data corre- 
sponding to 200 SPI exposures (^400 ks), we need to run 
3800 simulations (200 x 19). Using 32 10-core processors, 
the computing time is reduced from 950 to 3 days. 

4.3. Simulated Data Preparation 

The original data produced by Geant4 is a list of en- 
ergy deposits within sensitive volumes. These raw data 
require some selection and preparation before analysis. A 
low energy threshold of 20 keV is applied and the events 
occurring in dead detectors are removed. These reduc- 
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Figure 4. Comparison of the SPI spectral response and GEANT4 
simulation (blue) for a mono-energetic source (661 keV), after mass 
model improvements, for single events summed over all the detec- 
tors. The two curves, normalized to the photo-peak, differ by only 



tions affect the multiple events: some triple events will 
become doubles and some doubles will become singles. 
For the data analysis, we select the double events with 
energy deposits occurring in adjacent detectors. There 
are 42 corresponding pseudo detectors for all the possible 
couples of adjacent detectors. Then the simulated data 
are energy binned and recorded in the same format as 
the SPI data to simplify the analysis procedure. 
For a single pointing, simulations are carried out for 18 
P.A.s (100% polarized) and 1 unpolarized. However, the 
analysis requires simulated data for different polarization 
percentages for each angle. These data are produced by 
mixing the polarized simulated distributions with the un- 
polarized one using the formula 



G4(il» = 



IT x G40) (100 - n) x G4(f7) 



100 



100 



(2) 



where G4(LI, <p) represents the (Geant4) simulated data 
for a II percent polarized flux at the angle <f), G4(0) the 
100% polarized simulated data at angle <fr, and G4(C/) 
the unpolarized simulated data. 

4.4. Comparison between Geant^ Model and Standard 
SPI Response 

The search for polarization in a source emission re- 
quires accurate simulations of the instrument response, 
any differences would lead to a spurious polarization 
signal. Figure @] shows the spectral response of the 
GEANT4 sim ulation compared to the SPI standard re- 
sponse matrix ()Sturner et al.|[2003l ). for a mono-energetic 
source at 661 keV. The single events are summed over all 
the detectors and the counts per 1 keV bin are plotted in 
fraction of the photo-peak. With the implemented mass 
model improvements, the two curves differ by only ~ 8% 
in counts per bin with Compton to total counts ratio of 
0.32 and 0.35 for the GEANT4 simulation and the SPI 
response matrix, respectively. This ratio is important 
as it drives the production of multiples events. If the 
simulated detector geometry is too large or too heavy, 
this ratio would be smaller, i.e., less multiple events and 
more single events will be produced. Figure [5] shows the 
spatial distribution produced by the GEANT4 code com- 
pared to the measured SPI response, for the same config- 
uration. Compared to the SPI response, the unpolarized 
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Figure 5. Comparison of the SPI spatial distribution and 
GEANT4 simulation (blue and red) for the same exposure. The 
unpolarized simulation (blue) differs by only ~ 6% and the 20° 
polarized simulation (red) differs by ~ 20%. 

theoretical curve (blue) differs by only ~ 6% in counts 
per pseudo-detector whereas the 20° polarized one (red) 
differs by ~ 20%. It means that the signature of a po- 
larized signal is contained in the difference, i.e., ~ 15% 
of the recorded counts. The spatial distribution is highly 
modulated by the mask pattern and the telescope point- 
ing accuracy. Less obvious, the anti-coincidence system 
and especially its energy threshold has an important im- 
pact on the spatial distribution. A too high threshold 
will miss multiple events, increasing the recorded double 
events in the outer pseudo-detectors. 

5. DATA ANALYSIS 

5.1. Problem Formulation 

In a single SPI observation, the recorded flux is barely 
enough to produce an image or a spectrum of the source. 
This means that for a more complicated analysis such as 
polarization measurement, many observations need to be 
summed together. Unfortunately, the instrument con- 
figuration changes along the time: the pointing axis is 
shifted by 2°, every 30-40 minutes, while, on a longer 
timescale, detector failure must be accounted for together 
with other parameter (i.e., background) evolution. For a 
given source, our study relies on simultaneous analysis of 
a great number of pointings. Considering one exposure, 
the recorded data are modeled as 

D sd = xx G4 sd (n, 4>) + y x B sd , (3) 

where x is the source normalization, G4, the simulated 
count distribution, y is the background normalization 
and £>,the background spatial distribution. The source 
model G4 sc i (II, 4>) describes the number of counts for a 
science window s, in the pseudo detector d as a function 
of source polarization fraction, II, and angle, <j). The 
background spatial distribution is taken from an empty 
field observation temporally close to the observation. To 
be comparable with the data, the counts issued from the 
simulation are renormalized to the corresponding detec- 
tor livetimes. The x and y values are determined through 
a linear least squares resolution and the resulting \ 2 
value is stored. If the source and background flux can 
be considered constant on a long time scale, Equation 
([3]) is solved for the corresponding set of science win- 
dows. 

Note that it is different from previous polarization stud- 
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Figure 6. x 2 rnap produced during the data analysis. There are 19 X 201 x 2 values corresponding to the best fits over the P.A.s and 
P.F.s. A contour plot is projected above showing the confidence regions (68.27%, 95.45%, and 99.73%) around the best \ 2 ■ 



ies made with the SPI instrument, where the source and 
the background were not adjusted by a fitting procedure. 
For GRBs, the background was estimated using appro- 
priate time intervals before and after the source signal. 
This i s not feasib l e with most of the sources. For the 
Crab, iDean et all (|2008l ) considered a background with 
a constant pattern and its norm was calculated by sub- 
tracting the simulated counts from the data. This is 
hazardous because it implies the perfect knowledge of 
the instrument sensitivity and the source strength during 
each science window. In our code, the background pat- 
tern is built every six months from relevant empty field 
observations and its normalization is allowed to vary on a 
timescale from one science window to a revolution. Esti- 
mating the source and the background fluxes in the anal- 
ysis procedure introduces additional degrees of freedom 
(dof) but provides a better determination of the source 
contribution and allows polarization measurements for 
any (strong enough) source. 



5.2. Analysis Process 

The analysis process relies on the comparison between 
the recorded data and the theoretical ones. It can be 
described by the following procedure. 

• Science windows with the source closer than ±13° 
from the pointing direction are selected. This en- 
sures a good illumination of the detectors. 

• A quick standard analysis is run to filter the science 
window list from obvious data anomalies. 

• The SPI data are read into memory and the counts 
are summed over the selected energy range. 

• The corresponding simulated data are read into 
memory, the counts are summed over the selected 
energy range and weighted by the appropriate de- 
tector livetimes. 
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• If the source and the background fluxes can be 
considered as constant, the science windows are 
grouped per data sets (with a common x and y 
parameters). 

• Looping over the entire range of P.A.s <f> (0°-170°) 
and polarization fractions (P.F.s) II (—100% to 
100%), Equation (J3)) is solved simultaneously for 
each data set through a linear least-squares resolu- 
tion. 

• The corresponding x 2 ar e calculated, resulting in a 
19 x 201 x 2 ma P m the (0; n) space (see Figure[6]). 

• Using a cubic interpolation algorithm, the x 2 map 
is extended to 180 x 201. 

• The minimum x 2 i s used as the indicator of the 
best parameters. 

• A contour plot is displayed to visualize the con- 
fidence regions at the ler, 2a, and 3<r levels (see 
Figure EJ. 

Some additional information are extracted from the % 2 
map. 

Figure [7] shows the P.F. giving the best fit for each 
simulated angle (bottom) and the associated \ 2 (top). 

The best-fit solution appears in these curves at the 
lower x 2 value and at the maximum P.F., respectively. 
Note that only the P.F.s from 0% to 100% are physically 
correct. However, the analysis is run over negative po- 
larization percentages to find the best fits from a statis- 
tical point of view. Regarding the angular distribution 
predicted for polarized emissions (from the differential 
Compton cross-section), it is expected to have two math- 
ematical solutions (local minima): the first (good) one, 
with a positive fraction and the second, 90° away, with 
a negative fraction (understandable in terms of a pat- 
tern to be subtracted from a overestimated mean flux, 
see Figure |6]) . This symmetry is also seen in the polar- 
ization modulation (see Figure [7|) . 

5.3. Error Estimation 

The errors can be estimated from the A% 2 method 
depending on the n umber of parameters of interest 
(jLampton et al.lfl976| ). The la, 2a, and 3cr confidence 
interval for one parameter of interest, without regard to 
the value of the second, is given by Xmin + 1) Xmin + ^, 
and Xmin + 9- The errors on the P.A. and P.F. inde- 
pendently are calculated with these Ax 2 values. The 
confidence region (two dimensions), plotted above the 
X 2 map, is given by the x 2 criteria of two parameters 
of interest. The 68.27%, 95.45%, and 99.73% confidence 
regions are drawn at the level of Xmin + 2.30, Xmin + 6-18, 
and xLi n + 11-8. 

To be confident in our error analysis, we have verified 
this method by simulation using the Bootstrap method. 
It consists to run the same analysis thousands of times 
to quantify the errors. We apply the technique where the 
simulated counts are randomized within their (Gaussian) 
errors before to be re-analyzed. The best fit values of 
the P.A. and P.F. are stored for each trial. At the end, 
the pairs angle-fraction are spread around the real solu- 
tion and form a normal distribution (in two dimensions). 
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Figure 7. Results extracted from the x 2 rrmp. The P.F. giving 
the best fit for each position angle (bottom) and the corresponding 
X 2 (top). The errors are given at la 



The analysis of this distribution provides the region (pair 
angle- fraction) or the interval (angle or fraction) contain- 
ing 68.27%, 95.45%, and 99.73% of the realizations. We 
then identify on the x 2 map the A% 2 which gives the 
same confidence levels. We find that the 68.27%, 95.45%, 
and 99.73% confidence intervals are obtained at the level 
01 Xmin + 1> Xmin + 4 > and X m in + 9 > and the equivalent 
confidence regions at the level of x m i n + 2.30, x^ + 6.18, 
and Xmin + 11-8 as expected. 

6. APPLICATION 

6.1. The Crab Analysis 

The obvious candidate to search for polarization is the 
Crab (Pulsar and Nebula total emission) since it is one 
of the brightest source in the hard X-ray domain and 
the pulsar emission is expected to be highly polarized. 
Beyond that, a lot of polarization measurements exist 
at different wavelengths, giving a good set of results to 
compare with. 

For this analysis, we used 239 SPI science windows from 
revolutions 43 to 45 (2003 February), where the Crab 
Nebula was the main target. The photons fired in the 
simulations follow a power law distribution with a spec- 
tral index of 2.2 ranging from 100 keV to 3 MeV. We 
have considered the detected counts in the energy range 
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offering the best signal-to-noise ratio, 130 to 440 keV. 
The low energy threshold is fixed to ensure that Comp- 
ton efficiency is known with enough precision. The high 
energy one is related to the source signal-to-noise ratio 
decrease. 

The source and the background fluxes are considered as 
constant over one revolution (—3 days). The analysis 
process described in Section 15.21 is applied, producing a 
X 2 map in the polarization angle-fraction space. Figures 
[6]and[7l previously quoted for illustration purposes, come 
from this Crab observation analysis. The corresponding 
contour plot with the confidence regions is highlighted in 
Figure H The lowest x 2 (x 2 c d = im4 with dof= 10031) 
indicates a P.A. of 117° ± 9°, measured from north to 
east, with a P.F. of 28% ± 6% (errors given at la). 
The P.A., which indicates the direction of the pho- 
ton electric vector, appears to align nicely with the 
spin axis of the Crab pulsar, estimated at 124° ±0.1° 
from X-rays imaging (|Ng fc Rom ani 2004). Extending 
the energy range to 130-650 keV, we find a P.F. of 
32% ± 7% and 34% ± 8% in the energy range 130-1000 
keV. We did not find any significant dependence of the 
P.F. with the energy range. This behavior is consistent 
with a synchrotron emission process, which is expected 
for the Crab emission. Previous polarization measure- 
ments were made on the Crab, using only the off-pulse 
emission which is exp ected to be much more polarized. 
IKanbach etail (120051) repor t ed a P A. of 123° in the opti- 
cal domain and lForot etail (pOOl ) a P.A. of 122° ± 7.7° 
in the hard X-ray domain. Polarization measurements 
of the unpulsed emission has been already made with 
the SPI instrument, based on a previous mass model 
versio n and a different analysis process by iDean et al.l 
(2008). With -600 science windows (—1.5 Ms), covering 
the first three years of INTEGRAL operations (from rev- 
olution 43 to 422 in 2006 March), they found a P.A. of 
123° ± 11° with a P.F. of 46% ± 10%. Our work includes 
only 239 science windows (570 ks) and consider the total 
(pulsar and Nebula) emission, but we find a better con- 
strained solution, showing that the mass model improve- 
ments and the analysis process have strongly reduced the 
systematics errors. 

6.2. Systematic Errors 

To produce a spurious polarization signal requires 
an effect able to favor one direction in the scattered 
events, i.e., affecting more counts in some specific pseudo- 
detectors. This kind of systematic error could come from 
the instrument itself or/and from the simulations. Po- 
larimeters instruments are often designed to spin their 
detection plane during an observation so as to blend the 
potential detection artifacts. SPI is not spinning but its 
coded mask imaging capabilities require a specific ob- 
servation strategy (dithering) during which the pointing 
direction shift around the source. This means that the 
detectors, so the pseudo-detectors, lit by the source are 
not always the same. Moreover, when using many obser- 
vations over a one year period, the detection plane az- 
imuthal angle (rotation around the pointing axis) varies 
by 1 deg per day (mean value). These two characteristics 
drastically lim i t any instrumental systematic effect. 
iKalemci et al.l (|2007l ) have investigated the potential in- 
trinsic biases related to the detector with ground calibra- 
tion data. They could constrain them to be around 1%, 
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Figure 8. Contour plot from the analysis of the Crab Pulsar and 
Nebula. The best fit (x^ od = 1.034) indicates a P.A. of 117° ± 
9° with a P.F. of 28% ±6%. The 68.27%, 95.45%, and 99.73% 
confidence regions are shaded from dark to light gray. 

quite negligible for our study. 

On the other side, we have to check that the simula- 
tions are free of any systematic effects. To estimate these 
systematics, we have applied our analysis procedure to 
the long distance imaging ground ca libration data ob - 
tained at Bruyeres-le-Chatel, France (jAttie et al.l l2003"). 
We have selected 10 exposures (from 12 to 38 minutes) 
performed with an unpolarized source of 137 Cs at 125 
m. We have produced the corresponding simulated ob- 
servations with our model and selected the double events 
for which the total energy corresponds to the 137 Cs line 
energy. Applying our analysis process we found a P.F. 
< 5% whatever the angle. This value is overestimated 
due to the uncertainties linked to the calibration setup 
such as the position and the finite distance of the radioac- 
tive source. As a sec ond way, we have us ed the standard 
SPI response matrix (jSturner et a l. 2003) which provides 
the response for a given source position at infinity. From 
this standard response matrix, we have reproduced a set 
of simulated observations corresponding to the SPI ob- 
servations used in our Crab analysis. As this standard re- 
sponse matrix does not include the effects of polarization 
in the Compton scattering distribution, this simulated 
set is representative of observations of a non-polarized 
source. We have analyzed these observations with the 
procedure described above and found a P.F. < 4% what- 
ever the angle. Consequently, this demonstrates that the 
imaging system of SPI (pixelated detector + coded mask 
+ dithering) and the use of a complete and validated in- 
strument response ensure reliable polarization measure- 
ments with systematic effects smaller than 4%-5%. 

7. SUMMARY AND CONCLUSIONS 

We have presented specific analysis tools developed 
to study the polarization of emissions detected by the 
SPI spectrometer, aboard the INTEGRAL mission, 
in the hard X-ray domain (130 keV to 8 MeV for 
the polarimetry studies). Polarization information 
being contained in the Compton interaction parameters 
(diffusion and azimuthal angles), we used the photons 
which undergo a Compton interaction in one of the 19 
germanium crystals with the diffused photons escaping 
to an adjacent crystal. 

Precise simulations using Geant4 can provide the SPI 
instrument response for these polarized events, for 
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various P.A.s (from 0° to 180°) and P.F.s (from 0% to 
100%). The observed data are compared with each of 
these simulations, through \ 2 estimation in the same 
way than imaging is performed with SPI data. We 
produce a \ 2 map in the two-dimensional polarization 
angle-fraction space and consider the lowest \ 2 value to 
determine the most probable parameters of the incident 
flux. 

We have first applied our method to ^600 ks of obser- 
vation on the Crab and found immediately a P.A. of 
11 7° ±9°, in perfect agreement with previously publishe d 
value on the Crab (|Dean et alj|2008t iForot et al.ll2008f) . 
However, our ana lysis method is m ore precise than 
that presented in iDean et al.l (j2008j ) since we use an 
improved Geant4 version of the SPI/ INTEGRAL mass 
model and develop more sophisticated algorithms for the 
background determination and source flux extraction. 
In particular, we have included the SPI alignment 
correction with the spacecraft axis to improve the source 
projection into the detector plane, a parameter essential 
for the simulations. Moreover, our analysis process 
adjusts the source and background contributions within 
the fitting procedure and permits their variation in 
time. This allows us to better determine the source 
contribution and to apply the polarimetry analysis to 
other sources than strong GRBs or constant sources. 
In conclusion, the measurement of the polarization 
appears as a promising area in the hard X-ray domain, 
through the Compton interaction process inside de- 
tection systems. It requires a high statistic and thus 
strong sources. Nonetheless, pulsars, bright transients 
and GRB are as many potential targets for this kind 
of studies, which bring crucial information, very com- 
plementary to the spectral and timing analyzes. This 
promising and still unexplored field should provide key 
information on mechanisms at work inside the most 
powerful engines of the sky. 
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